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A few linkage maps of tea have been constructed using pseudo-testcross theory based on dominant marker 
systems. However, dominant markers are not suitable as landmark markers across a wide range of materials. 
Therefore, we developed co-dominant SSR markers from genomic DNA and ESTs and constructed a refer- 
ence map using these co-dominant markers as landmarks. A population of 54 Fi clones derived from recip- 
rocal crosses between 'Sayamakaori' and 'Kana-Ckl7' was used for the linkage analysis. Maps of both 
parents were constructed from the F t population that was taken for BCi population. The order of most of the 
dominant markers in the parental maps was consistent. We constructed a core map by merging the linkage 
data for markers that detected polymorphisms in both parents. The core map contains 15 linkage groups, 
which corresponds to the basic chromosome number of tea. The total length of the core map is 1218 cM. 
Here, we present the reference map as a central core map sandwiched between the parental maps for each 
linkage group; the combined maps contain 441 SSRs, 7 CAPS, 2 STS and 674 RAPDs. This newly con- 
structed linkage map can be used as a basic reference linkage map of tea. 
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Introduction 

Tea {Camellia sinensis) is one of the most popular beverage 
crops originated in Asia, and the products are widely con- 
sumed all over the world. In 2009, the FAO estimated that 
2.8 million hectares of tea fields was harvested (http:// 
faostat.fao.org/). Harvested fresh tea leaves are processed 
into many types of tea, including black tea, pan-fired green 
tea, steamed green tea, oolong tea and powdered green tea, 
and the processed tea leaves are traded around the world. 
Some of the major tea-producing countries are developing 
countries, for which tea is a major source of foreign currency 
revenue. Tea has a genome size of 4 Gb (Tanaka et al. 2006); 
its basic chromosome number is «=15 (Morinaga et al. 
1929) and it is commonly diploid. Since tea is an outcrossing 
woody plant with a self-incompatible reproductive system 
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(Tomo et al. 1956), its genome is always heterogeneous. 
Traditionally, tea trees have been propagated as seedlings, 
but seedling tea trees do not have synchronous bud break, 
uniform harvest times and uniform quality. Recently, clonal 
cultivars have become the main plant materials used in many 
countries. Clonal cultivars have synchronous bud break, uni- 
form harvest times and uniform quality, making them well 
suited for machine harvesting. 

To breed clonal tea cultivars, breeders have to select a 
large number of individuals with desirable characters from a 
large number of F] segregating populations. The selected in- 
dividuals are clonally propagated and screened with the in- 
formation of clonal strains. After evaluation, elite clones are 
released as cultivars. Production of high-quality tea leaves is 
the most important breeding objective, but the evaluation is 
laborious and time-consuming. In addition to being a 
lengthy process, tea breeding also requires large field areas. 
For these reasons, it is efficient and effective to multiply 
early-stage populations and select breeding materials on 
the basis of an evaluation of DNA markers for important 
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agronomic characters that are controlled by a few loci, 
such as pest resistance, before planting the seedlings in the 
field. At Makurazaki Tea Research Station, a breeding sys- 
tem that combines marker-assisted selection (MAS) in the 
juvenile phase and the acceleration of generations has shown 
successful results (Tanaka 2006). It would be efficient to be 
able to perform MAS for many characters in many different 
populations, so a method of genetic analysis is needed for a 
wide range of genetic materials. It is therefore desirable to 
construct a reference map with many landmark markers. 

A few linkage maps of tea based on the pseudo-testcross 
theory (Grattapaglia and Sederoff 1994) have been reported 
(Hackett et al. 2000, Ota and Tanaka 1999, Tanaka et al. 
1995). These maps were constructed by using dominant 
marker systems such as Random Amplified Polymorphic 
DNAs (RAPDs) or Amplified Fragment Length Polymor- 
phisms (AFLPs). By using pseudo-testcross theory, it is easy 
to construct accurate linkage maps for small F] populations 
having two heterozygous parents, because the segregation 
ratio of each marker is 1 : 1 and linkage relationships among 
them are very simple and clear. However, the utility of dom- 
inant markers depends on the particular material being 
tested: they are not universal markers. In addition, the rela- 
tionships among the linkage groups on the previously report- 
ed maps are unknown. In a few fruit trees, reference linkage 
maps have been developed by using co-dominant Cleaved 
Amplified Polymorphic Sequences (CAPS) and Simple 
Sequence Repeats (SSR) markers (Omura et al. 2000, 
Yamamoto et al. 2002, 2005, 2007). Here, we developed a 
high-density reference map of tea, using both newly devel- 
oped SSR markers and some previously reported SSR and 
CAPS markers (Ogino et al. 2009, Taniguchi et al. 2012, 
Ujihara et al. 2011) as landmark markers. RAPD markers 
were included in the parental maps but were not used as 
landmark markers because of their dominant nature. 

Materials and Methods 

Plant materials 

An F, segregating population consisting of 54 clones was 
derived from reciprocal crosses between the Japanese green 
tea cultivar 'Sayamakaori' and the strain 'Kana-Ckl7'. 
'Sayamakaori' was selected from a population derived from 
naturally pollinated seedlings of 'Yabukita' (Fuchinoue et 
al. 1972); the pollen parent of 'Sayamakaori' is believed to 
be germplasm introduced from China (Tanaka et al. 2001). 
'Kana-Ckl7' is germplasm introduced from China. It has 
been reported that Chinese tea materials have greater genetic 
diversity than Japanese ones (Kaundun and Matsumoto 
2003, Wachira et al. 2001, Yamaguchi and Tanaka 1999). 
We expected there to be high heterozygosity across the ge- 
nomes of these materials because of the lack of duplication 
in the pedigrees of the parents and because tea has an out- 
crossing propagation system. 



Development of SSR markers 

The SSR markers with names beginning with 'TMSLA' 
and 'MSG' (Supplemental Table 1) were developed from 
microsatellite-enriched libraries. These libraries, enriched 
for GA and GT, were developed by Genetic Identification 
Services Inc. (GIS, Chatsworth, CA, USA) from DNA of 
strain 'KM62', extracted by using a modified CTAB method 
(Tanaka et al. 2001). Sequencing was performed by using 
a model 310, 3130x/, or 3700 Genetic Analyzer (Applied 
Biosystems, Branchburg, NJ, USA). The primer sets were 
designed with Primer Express (Applied Biosystems) or 
read2Marker (Fukuoka et al. 2005) software. The SSR 
markers with names beginning with 'MSE' (Supplemental 
Table 1) were derived from EST databases; the correspond- 
ing primer sets were also designed with read2Marker. SSR 
marker TMS3 was developed by using a PCR method 
(Tanaka and Ikeda 2000). 

SSR marker detection 

DNA was extracted from fresh leaves by a method using 
diatomaceous earth and spin filter (Tanaka and Ikeda 2002). 
The PCR conditions were similar to those previously de- 
scribed (Tanaka et al. 2003). PCR reactions (10 ul total 
volume) contained 0.05 units/ ul of AmpliTaq DNA Poly- 
merase (Applied Biosystems), the attached reaction buffer, 
2.0 mM MgCl 2 , 0.16mM of each dNTP, 1% formamide, 
2 ng/ul template DNA and 2 uM of each primer DNA. The 
PCR temperature conditions were based on the touchdown 
PCR technique (Don et al. 1991). The PCR program was as 
follows: 5 minutes at 94°C to completely denature the DNA; 
followed by 34 cycles of 30 s at 94°C, 60 s at annealing tem- 
perature (described below) and 30 s at 72°C; followed by 
10 minutes at 72°C to allow complete double-strand DNA 
synthesis. The annealing temperature, which was 62°C in 
the first cycle, was reduced 0.5°C per cycle during cycles 2 
to 14 and maintained at 55°C for the last 20 cycles. PCR 
products were post-PCR labeled (Kukita and Hayashi 2002), 
then separated and detected by using a model 3 1 0, 3 1 3 0x1, or 
3700 Genetic Analyzer (Applied Biosystems). The sizes of 
the amplified fragments were calculated on the basis of 
internal-standard DNA (Gene Scan 500 LIZ, Applied Biosys- 
tems) by using GeneMapper software (Applied Biosystems). 

CAPS marker detection 

Total DNA was extracted by using a modified CTAB 
method (Kaundun and Park 2002). PCR, restriction diges- 
tion and electrophoresis were performed according to the 
methods described by Ujihara et al. (2005). CAPS primer 
sequences are given in Supplemental Table 2. 

RAPD marker detection 

Template DNA for RAPD marker detection was prepared 
in the same way as for SSR analysis. The PCR reactions 
(10 ul total volume) contained 0.5 units of Ampli-Taq DNA 
Polymerase (Applied Biosystems), the attached reaction 
buffer, 2.5 mM MgCl 2 , 0.125 mM of each dNTP, 2 ng/ul 
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template DNA and 2 uM primer DNA. The PCR program 
was as follows: 5 minutes at 93°C to completely denature the 
DNA; followed by 40 cycles of 1 min at 93 °C for denatura- 
tion, 1.5 min at 42°C for annealing and 1 min at 72°C for 
extension; followed by 10 min at 72°C to allow complete 
double-strand DNA synthesis. The PCR products were sep- 
arated by electrophoresis in 2% agarose gel and visualized 
by using a UV transilluminator after being stained with 
ethidium bromide. The name of a RAPD marker refers to the 
primer name and expected band size. For example, 
'OpA_12_600' represents the marker band of 600 bp ob- 
tained with primer OpA-12 ('Op' means Operon Technolo- 
gies, Inc. Alameda, CA, USA). 

Mapping algorithm 

The F] segregating population we used in this analysis 
was not very large, so we formulated and applied a plan to 
develop a three-line reference map by using information 
from both parents along with the F] (Fig. 1). This reference 
map was designed to provide an accurate marker orders and 
linkage distances between each markers, even by with limit- 
ed population size. AntMap (Iwata and Ninomiya 2006) 
software, which is based on ant-colony optimization theory 
(Dorigo et al. 1996), was used to construct the linkage maps. 
On the basis of the double pseudo-testcross theory 
(Grattapaglia and Sederoff 1994), the F| segregating popula- 
tion was regarded as resembling a BCi population for the 
purpose of linkage analysis. Before the linkage analysis, the 
segregation data were sorted into data sets derived from each 
parent (Fig. 1). Then, a map of each parent was constructed 
by using the sorted data sets. In many cases, SSR marker 
bands generated by the same primer sets were expected to 
amplify the same locus in each of the two parents, so SSR 
markers tended to map within linkage groups in the same 
order in both parental maps. Thus, the amount of genotype 
data was doubled for those markers that detected polymor- 
phisms within both parents. We manually selected those 
common markers, which were assumed to be derivatives 
from the same locus in both parents because of the segrega- 
tion patterns of the flanking markers. This selected marker 
set was designated as the core marker set. We then analyzed 
the linkage relationships among the only core marker set on 
the basis of the combined (54 + 54= 108) genotyping data 
and a new map was constructed that consisted of only the 
core marker set. This merged map was designated as the 
core map. To present the relationship between the linkage 
groups of the core map and the two parental maps, we illus- 
trated the reference map in such a way that the core map for 
each linkage group was sandwiched between the parental 
maps by using MapChart software (Voorrips 2002). 

RAPD STS 

Methylene-blue-stained agarose gel blocks containing 
target RAPD bands were cut out and re-amplified by PCR 
for 1 5 cycles as described above for RAPD marker detec- 
tion. Re-amplified PCR products were ligated into the 



pGEM-T vector (Promega, Madison, WI, USA). The plas- 
mid was used to transform E. coli strain DH5a and se- 
quenced with a 3130x/ Genetic Analyzer (Applied Biosys- 
tems). 

To clone the flanking sequences of RAPD STS markers, 
thermal asymmetric interlaced polymerase chain reaction 
(TAIL-PCR, Liu et al. 1995) was performed using arbitrary 
primers MAT01 (CNSCTSCTNTWTCTT), MAT02 (CNNC 
WNCAANAWCAA), MAT03 (GNSGASGANAWAGAA), 
MAT04 (GNNGWNGTTNTWGTT) and specific primers 
designed by using STS sequence information. The PCR so- 
lutions (10 (al total volume) were the same as for SSR mark- 
er detection except for the arbitrary primer concentration 
(50 uM). The temperature program basically followed the 
TAIL-PCR protocol (Liu et al. 1995), integrating the touch- 
down PCR technique. The annealing temperature for the first 
PCR, which was 60°C during the first cycle, was reduced 
1°C per cycle during the first 5 cycles. In the asymmetric 
super-cycle, the temperature for non-specific annealing was 
maintained at 38°C and the temperature for specific anneal- 
ing started at 60°C and was reduced by 0.4°C per cycle for 
15 super-cycles. For the second and third PCRs, the temper- 
ature for non-specific annealing was maintained at 38°C and 
the temperature for specific annealing started at 60°C and 
was reduced by 0.5°C per cycle for 12 super-cycles. For the 
third PCR, primers with single nucleotides (A, T, G or C) 
added to the 3 '-end of the arbitrary primer sequences were 
used for selective amplification of the flanking sequence. 

Results 

Development of SSR markers 

A total of 733 primer sets were designed from 1440 
clones from SSR-enriched libraries and 298 primer sets were 
designed from 5504 tea unigene sequences obtained in our 
laboratory and from publicly available EST information. A 
total of 377 primer sets generated polymorphic amplicons 
between two alleles within one or both parents. The newly 
developed SSR markers are described in Supplemental 
Table 1. 

Construction and characterization of linkage maps 

Three-line reference linkage maps of the tea clonal cul- 
tivar 'Sayamakaori' and clonal strain 'Kana-Ckl7' were 
constructed using SSR, RAPD and CAPS markers (Fig. 2). 
Details of the markers on the maps are described in 
Supplemental Tables 1, 2. The map of 'Sayamakaori' con- 
sisted of 701 loci, 363 of which are SSR markers. Similarly, 
the map of 'Kana-Ckl7' consisted of 701 loci, 354 of which 
are SSR markers. Both maps contained 15 linkage groups, 
which are summarized in Table 1 . All 441 SSR markers with 
normal segregation within the Fi population could be 
mapped to one of the 15 linkage groups. The maps of 
'Sayamakaori' and 'Kana-Ckl7' encompassed a genetic dis- 
tance of 1305 cM and 1298 cM, respectively, with an aver- 
age distance between markers of 1.93 cM (Sayamakaori) 
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DNA marker segregation patterns on putative one linkage group in a segregating population derived from crossing Parent A and Parent B 
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Fig. 1. Pattern diagram of the construction of three-line reference map based on pseudo-testcross theory. 
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Fig. 2. Three-line reference genetic linkage map of tea (Camellia sinensis) using an F[ population derived from reciprocal crosses between 
'Sayamakaori' (SAK) and 'Kana-Ckl7' (KnCk). Locus names are given to the right of each map; genetic distances (in cM) are to the left. Core 
markers are denoted by bold characters. 
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Fig. 2. (continued) 



and 1.85 cM (Kana-Ckl7). These two maps each included 
the 279 SSR markers and 2 CAPS markers in the core mark- 
er set. 

The preliminary core map consisting of SSR and CAPS 
core markers had a large gap region in linkage group 5. To 
supply a marker to cover this gap, RAPD marker 
'OpT_18_200' f rom the map of 'Kana-Ckl7' was cloned 
and sequenced. Because there is no polymorphism for this 
sequence in 'Sayamakaori', the flanking sequence was am- 



plified by using TAIL-PCR and sequenced. There were 2 
SNPs revealing heterozygosity in 'Sayamakaori' in the 
flanking sequence of 'OpT_18_200' (accession numbers: 
AB623060 [Sayamakaori], AB623061 [Kana-Ckl7]); these 
SNPs were mapped to the same locus as RAPD marker 
'OpJVl l 1 750' on the 'Sayamakaori' map. The segregation 
patterns obtained for 'OpM_l 1_750' and the new marker 
were identical, so we mapped the new STS marker 
(designated 'OpM_18_200STS') onto the core map. 
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The revised core map contained 15 linkage groups with a 
total length of 1218 cM. The average distance between 
markers was 4.35 cM, without any gap region. 

Discussion 

Characterization of the newly constructed linkage maps 

The core map and the parental maps each covered 15 
linkage groups; this corresponds to the basic chromosome 
number (n= 15) of tea. All 441 SSR markers with normal 
segregation within this population could be mapped to one 
of the 15 linkage groups. The total length of the core map is 
1218 cM, and that of the combined (three-line) map is 
1317 cM. The genome size of tea is estimated to be 4 Gbp 



(Tanaka et al. 2006), so 10 cM of genetic distance would 
correspond to 30 Mbp of physical distance on average. 

Illustration and application of the reference linkage map 

We have illustrated the reference map in such a way that 
the core map for each linkage group is sandwiched between 
the corresponding parental maps. Although obtained from 
the same Fi population, the sets of linkage information for 
each parent are independent. If the maps of the two parents 
were to be merged into a single map by using the landmark 
markers, the order of the non-landmark markers would be 
expected to be very inaccurate. This three-line manner of 
illustration presents reasonably accurate marker orders. 
It is common to have inaccurate marker orders in linkage 
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Table 1. Summarized data of each linkage group of tea (Camellia sinensis) reference map 







Genetic distance 
(cM) 


Total distance 
(cM)" 


Number of loci 


SSR loci 


RAPD loci 6 


CAPS and 
STS loci 


Linkage group 1 


Core 


110 


120 


18 


18 
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120 
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21 


22 






Kana-Ckl7 
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44 


27 


17 




Linkage group 2 


Core 
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106 


17 


17 








Sayamakaori 


102 




32 


17 


15 






Kana-Ckl7 


106 




40 


22 


18 




Linkage group 3 


Core 


97 


100 


23 


23 








Sayamakaori 


101 




62 


31 


28 


3 




Kana-Ckl7 


111 




60 


32 


28 




Linkage group 4 


Core 


98 


98 


21 


21 








Sayamakaori 


108 




52 


26 


25 


1 




Kana-Ckl7 


72 




55 


28 


27 




Linkage group 5 


Core 


67 


97 


13 


12 




1 




Sayamakaori 


94 




37 


15 


21 


1 




Kana-Ckl7 


86 




38 


16 


26 


1 


Linkage group 6 


Core 


86 


92 


31 


31 








Sayamakaori 


92 




85 


45 


40 






Kana-Ckl7 


88 




75 


34 


41 




Linkage group 7 


Core 

Sayamakaori 


84 
92 


91 


14 
24 


14 
16 


8 






Kana-Ckl7 


86 




31 


19 


12 




Linkage group 8 


Core 


90 


90 


31 


30 




1 




Sayamakaori 


96 




69 


41 


27 


1 




Kana-Ckl7 


97 




59 


39 


19 


1 


Linkage group 9 


Core 


85 


87 


23 


23 








Sayamakaori 


71 




52 


32 


19 


1 




Kana-Ckl7 


97 




60 


30 


30 




Linkage group 10 


Core 

Sayamakaori 


78 
86 


84 


16 
43 


16 
21 


22 






Kana-Ckl7 


90 




45 


19 


26 




Linkage group 1 1 


Core 


74 


79 


15 


15 








Sayamakaori 


77 




29 


19 


10 






Kana-Ckl7 


77 




35 


22 


13 




Linkage group 12 


Core 


60 


75 


16 


16 








Sayamakaori 


58 




39 


23 


16 






Kana-Ckl7 


80 




48 


21 


26 


1 


Linkage group 13 


Core 


62 


73 


15 


15 








Sayamakaori 


65 




47 


21 


26 






Kana-Ckl7 


70 




38 


16 


22 




Linkage group 14 


Core 


62 


64 


20 


20 








Sayamakaori 






oz 


TO 
ZO 








Kana-Ckl7 


58 




47 


22 


25 




Linkage group 1 5 


Core 


61 


61 


6 


5 




1 




Sayamakaori 


57 




25 


7 


17 


1 




Kana-Ckl7 


73 




26 


7 


18 


1 


Total 


Core 


1218 


1317 


279 


276 


0 


3 




Sayamakaori 


1305 




701 


363 


330 


8 




Kana-Ckl7 


1298 




701 


354 


348 


4 


Total marker number without duplication 






1123 


441 


678 


9 



° Total distance is culculated by Core linkage group plus extra region of parent maps. 
* The number of RAPD loci does not include RAPD-STS. 
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maps because of missing data, especially when the popula- 
tion size is small. In this study, we used automated DNA se- 
quencers to detect the sizes of amplified fragments from the 
SSR marker analysis. In this method, some background 
peaks hampered the detection of particular alleles; for exam- 
ple, this occurred with MSG0213 in linkage group 4 of 
'Sayamakaori' and MSE0230a in linkage group 3 of 'Kana- 
Ckl 7'. For such markers, the core map should be more accu- 
rate than either parental map because it is based on double 
the amount of linkage data per marker. 

In addition, the core markers are likely to be more infor- 
mative and more useful than other markers for genetic anal- 
ysis in other populations. For a marker to be designated as a 
core marker, both parental materials had to be heterozygous 
for the marker and the segregation of the alleles derived 
from one parent had to be detectable without the interference 
of alleles from the other parent. That is, if the genotype of 
one parent was AB, the genotype of the other parent had to 
have been AC, BC or CD. In other words, a core marker 
must have had at least three distinguishable alleles within 
the two parental lines. To select markers for genetic analy- 
sis, it is efficient to choose markers firstly from the core 
marker set, secondly co-dominant markers on the parental 
maps. If there are no polymorphic markers within the target 
area, it is possible to search SNP markers from the sequence 
information of SSR markers. Even if no SSR markers are 
found in the target region, it is possible to search SNPs from 
the sequence of flanking region of SSR markers, or even 
from the sequence of an RAPD band or its flanking region, 
as was done here in the case of 'OpT_18_200' STS. 

The linkage map presented here should be useful as a ref- 
erence map of tea. This reference map will play an important 
role in establishing the alignment of genomic sequences of 
tea as we work toward full-genome sequencing. In the fu- 
ture, we expect that many additional markers will be devel- 
oped and mapped and these will further advance the accura- 
cy and density of the available maps. 

Acknowledgments 

The authors thank Mrs. M. Iwata, Mrs. Y. Kuramae and 
Mrs. N. Rokutan, who are assistants at the Makurazaki Tea 
Research Station of NARO Institute of Vegetable and Tea 
Science, for their help. We also thank Dr. H. Kato of the 
NARO Institute of Crop Science and Dr. T. Hayashi of 
National Agricultural Research Center for their helpful com- 
ments on the manuscript. This work was supported by a 
grant from the Ministry of Agriculture, Forestry and Fisher- 
ies of Japan (Rice Genome Project DM-1503) and by NARO 
Research Project No. 211, 'Establishment of Integrated 
Basis for Development and Application of Advanced Tools 
for DNA Marker- Assisted Selection in Horticultural Crops.' 

Literature Cited 

Don.R.H, B.J.Cox, B.J.Wainwright and J.S.Mattick (1991) 'Touch- 



down' PCR to circumvent spurious priming during gene amplifica- 
tion. Nucleic Acids Res. 19: 4008. 

Dorigo,M., V.Maniezzo and A.Colorni (1996) The ant system: optimi- 
zation by a colony of cooperation agents. IEEE Trans. Syst. Man. 
Cybern. B26: 29-41. 

Fuchinoue,Y., M. Tanaka, S.Funakoshi and N.Okano (1972) A newly 
registered variety 'Sayamakaori' for green tea. Tea Res. J. 37: 12- 
23. 

Fukuoka,H., T.Nunome, Y.Minamiyama, I. Kono, N.Namiki and A. 
Kojima (2005) read2Marker: a data processing tool for microsatel- 
lite marker development from a large data set. BioTechniques 39: 
472^176. 

Grattapaglia, D. and R. Sederoff (1994) Genetic linkage maps of 
Eucalyptus grandis and Eucalyptus urophylla using a pseudo- 
testcross: Mapping strategy and RAPD markers. Genetics 137: 
1121-1137. 

HacketfC.A., F.N.Wachira, S.Paul, W.Powell and R.Waugh (2000) 
Construction of a genetic linkage map for tea (Camellia sinensis). 
Heredity 85: 346-355. 

Iwata,H. and S.Ninomiya (2006) AntMap: Constructing genetic link- 
age maps using an ant colony optimization algorithm. Breed. Sci. 
56: 371-377. 

Kaundun,S.S. and Y.G.Park (2002) Genetic structure of six Korean 
tea plantations as revealed by RAPD-PCR markers. Crop Sci. 42: 
594-601. 

Kaundun,S.S. and S.Matsumoto (2003) Development of CAPS mark- 
ers based on three key genes of the phenylpropanoid pathway in 
tea, Camellia sinensis (L.) O. Kuntze, and differentiation between 
assamica and sinensis varieties. Theor. Appl. Genet. 106: 375-383. 

Kukita,Y. and K. Hayashi (2002) Multicolor post-PCR labeling of 
DNA fragments with fluorescent ddNTPs. BioTechniques 33: 502- 
506. 

Liu,Y.G, N.Misukawa, T.Oosumi and R.F.Whittier (1995) Efficient 
isolation and mapping of Arabidopsis thaliana T-DNA insert junc- 
tions by thermal asymmetric interlaced PCR. Plant J. 8: 457^163. 

Morinaga,T., E.Fukushima, T.Kano, Y.Maruyama and Y.Yamasaki 
(1929) Chromosome numbers of cultivated plants II. Bot. Mag. 
Tokyo 43: 589-594. 

Ogino,A., J. Tanaka, F. Taniguchi, M.R.Yamamoto and K.Yamada 
(2009) Detection and characterization of caffeine-less tea plants 
originated from interspecific hybridization. Breed. Sci. 59: 277- 
283. 

Omura,M., T.Ueda, M.Kita, A.Komatsu, Y.Takanokura, T.Shimada, 
T.E.Inagaki, H.Nesumi and T.Yoshida (2000) EST Mapping of 
Citrus. Proceedings of the International Society of Citriculture IX 
Congress Orlando, Fla, USA: 1: 71-74. 

Ota,S. and J. Tanaka (1999) RAPD-based linkage mapping using F! 
segregating populations derived from crossings between tea culti- 
var 'Sayamakaori' and strain 'Kana-Ckl7'. Breed. Res. 1 (Suppl. 
1): 16. 

Tanaka, J. (2006) Breeding system combining of juvenile selection us- 
ing DNA markers and acceleration of generations, aiming for the 
significant improvement of the tea breeding in Makurazaki Tea 
Research Station, NIVTS. Breed. Res. 8: 1 19-126. 

Tanaka,J., Y.Sawai and S. Yamaguchi (1995) Genetic analysis of 
RAPD markers in tea. Jpn. J. Breed. 45 (Suppl. 2): 198. 

Tanaka, J. and S.Ikeda (2000) Generation of microsatellite markers by 
PCR in tea. Breed. Res. 2 (Suppl. 1): 37. 

Tanaka,J., N. Yamaguchi and Y.Nakamura (2001) Pollen parent of tea 
cultivar Sayamakaori with insect and cold resistance may not exist. 
Breed. Res. 3: 43^18. 



High-density reference linkage map of tea 



273 



Tanaka,J. and S.Ikeda (2002) Rapid and efficient DNA extraction 
method from various plant species using diatomaceous earth and a 
spin filter. Breed. Sci. 52: 151-155. 

Tanaka,J., S.M.Ota and Y.Takeda (2003) The garden-variety camellia 
'Robiraki' derived from crossing between Camellia japonica as a 
seed parent and C. sinensis as a pollen parent. Application of 
RAPD and an SSR marker analysis to tea breeding by interspecific 
hybridization. Breed. Res. 5: 149-154. 

Tanaka,J., F.Taniguchi, N.Hirai and S.Yamaguchi (2006) Estimation of 
the genome size of tea (Camellia sinensis), camellia (C. japonica), 
and their interspecific hybrids by flow cytometry. Tea Res. J. 101: 
1-7. 

Taniguchi,F., H.Fukuoka and J.Tanaka (2012) Expressed sequence 
tags from organ-specific cDNA libraries of tea (Camellia sinensis) 
and polymorphisms and transferability of EST-SSRs across 
Camellia species. Breed. Sci. 62: 186-195. 

Tomimoto,Y., H.Ikehashi, K.Kakeda and Y.Kowyama (1999) A pistil- 
specific PR-1 like protein of Camellia, its expression, sequence 
and genealogical position. Breed. Sci. 49: 97-104. 

Tomo,J., Y.Fuchinoue and H.Fuchinoue (1956) Study of self- 
fertilization of the tea plant. Tea Res. J. 7: 14-20. 

Ujihara,T., S.Matsumoto, N.Hayashi and K.Kohata (2005) Cultivar 
identification and analysis of the blended ratio of green tea produc- 
tion on the market using DNA markers. Food Sci. Technol. Res. 



11:43-45. 

Ujihara,T., F.Taniguchi, J.Tanaka and N.Hayashi (2011) Develop- 
ment of expressed sequence tag (EST)-based cleaved amplified 
polymorphic sequence (CAPS) markers of tea plant and their appli- 
cation for cultivar identification. J Agric. Food Chem. 59: 1557- 
1564. 

Voorrips,R.E. (2002) MapChart: software for the graphical presenta- 
tion of linkage maps and QTLs. Heredity 93: 77-78. 

Wachira,F.N., J.Tanaka and Y.Takeda (2001) Genetic variation and 
differentiation in tea germplasm revealed by RAPD and AFLP 
variation. J. Hortic. Sci. Biotech. 76: 557-563. 

Yamaguchi,N. and J.Tanaka (1999) Classification of Japanese original 
tea by maternal inherited RAPDs. Breed. Res. 1 (Suppl. 2): 295. 

Yamamoto,T., T.Kimura, M.Shoda, T.Imai, T.Saito, Y.Sawamura, 
K.Kotobuki, T.Hayashi and N.Matsuta (2002) Genetic linkage 
maps constructed by using an interspecific cross between Japanese 
and European pears. Theor. Appl. Genet. 106: 9-18. 

Yamamoto,T., M.Yamaguchi and T.Hayashi (2005) An integrated ge- 
netic linkage map of peach by SSR, STS, AFLP and RAPD. J. Jpn. 
Soc. Hort. Sci. 74: 204-213. 

Yamamoto,T., T.Kimura, S.Terakami, C.Nishitani, Y.Sawamura, T. 
Saito, K.Kotobuki and T.Hayashi (2007) Integrated reference ge- 
netic linkage map of pear based on SSR and AFLP markers. Breed. 
Sci. 57: 321-329. 



